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intermediate RHIC energy, ^/s^ = 62.4 GeV, are computed at large p in pertur- 

bative QCD. Next-to-leading order calculations in p-p scattering are first presented. 

The quenching of the prompt photon and pion yield due to energy loss processes 

in central Au-Au with respect to p-p collisions is then predicted. At this energy, 

the small phase-space available to produce hard partons makes the pion quenching 

almost as pronounced as at v /s NN = 200 GeV, despite the smaller gluon density of 

the produced medium. In the photon sector, energy loss effects prove small because 

of this very phase-space restriction, which favours the direct production channel. A 

significant suppression of high-p ± photons is however predicted, because of a strong 

isospin effect together with the depletion of nuclear parton densities at large x. 
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1. Introduction 

Undoubtedly, the strong quenching of leading hadrons up to p ± ~ 20 GeV in central 
Au-Au collisions at v /s NN = 200 GeV is one of the most significant results obtained 
so far at RHIC [1]. Remarkably, these data appear to be consistent with the "jet- 
quenching scenario" [2], where the high-p ± hadron depletion is due to energy loss 
processes coming from the parton multiple scattering in a dense QCD medium [3]. 
Despite these spectacular experimental results, however, our current theoretical un- 
derstanding of the fragmentation process modified in QCD media is still not clear. 
Consequently, there is an important need, on the phenomenological side, to investi- 
gate possible energy loss effects on a variety of observables. 

Among others, the production of prompt photons is a particularly interesting 
probe (see e.g. [4]). For quite some time, the usual belief had been that prompt 
photons, because of colour neutrality, should not be affected by the dense medium 
produced in high-energy nuclear collisions and, hence, could serve as a reliable base- 
line for "coloured" hard probes. This, however, may not be correct, and was first 
called into question a few years ago by Jalilian-Marian, Orginos and Sarcevic [5]. 
Indeed, in QCD perturbation theory at leading-order, prompt photons are produced 
in the partonic process, as in the Drell-Yan mechanism for instance, but also by the 
collinear fragmentation of hard quarks and gluons [6], just like any leading hadron. 
There is thus no reason, a priori, why the dense medium responsible for the significant 
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quenching of hadrons should not affect the latter channel as well. As a matter of fact, 
the PHENIX collaboration recently reported on a preliminary measurement of the 
prompt photon quenching factor in central Au-Au collisions at y/s NN = 200 GeV [7], 
which may be compatible with energy loss effects in the photon sector [8]. 

Another promising possibility to understand the dynamics of dense medium for- 
mation is to investigate, in a systematic way, the energy dependence of the quench- 
ing phenomenon [9]. In particular, the preliminary data in Au-Au collisions at 
y/s NN = 62.4 GeV, together with the ongoing analysis of p-p collisions at the same 
energy (Run-6), may actually provide a missing step between SPS [10] and RHIC 
results [1]. 

The aim of this paper is therefore to investigate the quenching of prompt photons 
(as well as that of hard pions for consistency) at intermediate RHIC energy, y/s NN = 
62.4 GeV. Predictions at this energy will be given in Section |3] and compared with 
estimates at top RHIC energy, ^/s NN = 200 GeV. Before discussing these results, we 
first detail the model used to compute photon and pion spectra in p-p and Au-Au 
collisions in the next Section. 



2. Model 

2.1 Perturbative QCD framework in p—p collisions 

Single-inclusive pion and photon hadroproduction cross sections in p-p collisions 
are computed in QCD at next-to- leading order (NLO) in a s , using the work of 
Refs. [6,11,12]: 

d d ^ d = J dxi dx2 F i/ P ( X ^ M ) F j/p (x 2 ,M) ^ D„ /k (z,M F ) 
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for photons. In the latter case, we marked the distinction between the photons pro- 
duced directly (labelled "dir" ) in the partonic scattering process, and those produced 
from the collinear fragmentation of the hard parton k (labelled "frag"). 

In Eqs. ( j2.1|) and (j2.2|), Fi/ P stands for the parton distribution functions (PDF) of 
flavour i in a proton, D n /k (D 1 /k) for the fragmentation functions (FF) of the parton 
k into a pion (photon), and (respectively, Ky^) for the leading-order (respec- 
tively, next-to- leading order) partonic cross section computed in the MS scheme 1 . We 
denote by /i, M and Mp the renormalisation, the factorisation and the fragmentation 
scale. 

In this study, we use the NLO CTEQ6M parton density in a proton [13], and the 
NLO Kniehl-Kramer-Potter (KKP) [14] 2 and the NLO Bourhis-Fontannaz-Guillet 
(BFG) [16] fragmentation functions into a pion and into a photon, respectively. 
Following the scale-fixing procedure used in [8], scales are chosen so as to minimise 
the scale-dependence of the NLO predictions. It turns out that rather low scales fulfil 
this requirement. We choose in the following the range p 2 ± /8 < [i 2 = M 2 = M 2 < p 2 ± , 
both for pions and photons, as a measure of the uncertainty of the NLO calculations. 

In the following, we shall be mostly interested in the nuclear effects on single 
pion and photon p x -spectra in the 5% most central Au-Au collisions. Therefore, 
only quenching factors 



are considered, where <rf e ° is the geometric cross section obtained via the Glauber 

' AuAu ° 

multiple scattering theory, <t nn the nucleon-nucleon cross section, and (A^ co n)| c the 
number of binary nucleon-nucleon collisions in central Au-Au collisions [4]. The 
ratio is properly normalised so that any "nuclear" effect (in a broad sense) in Au-Au 
scattering would make the quenching factor deviating from 1. Moreover, because 
of the present lack of understanding of medium-modified fragmentation processes, 
Au-Au spectra (and thus p-p spectra in Eq. (|2.3|) ) are determined at leading-order 
accuracy. 

2.2 Nuclear effects in Au— Au collisions 

Several nuclear effects need to be considered when computing pion and photon large- 
p ± spectra in Au-Au collisions. The most obvious one concerns the replacement of 



1 We omit the explicit dependence of &ij t k and -Ky.fc on the kinematic variables x\, Xi, y/s, p ± , 
and y for clarity. 

2 Albino, Kniehl, and Kramer (AKK) provide a more recent set of fragmentation functions [15]. 
We use, however, the KKP fit for consistency with Ref. [8]. 
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the proton PDF in Eqs. Q) and Q by: 



F i/A (x, M) = Z F i/P (x, M) + (A- Z) F i/n {x, M), 
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for a nucleus with Z protons and A — Z neutrons. The neutron PDF Fj/„ in Eq. (|2.4| ) 
is obtained from the proton one, Fi/ P , by isospin conjugation: u p = d n , d p = u n , 
yv = d n , dP = u n , s p = s n , and g p = g n . 

In addition to the above trivial isospin effect, nuclear parton densities differ from 
those probed in proton targets, over the whole x range (see [17] for a recent review). 
In order to take into account such a "genuine" nuclear effect (as opposed to the 
isospin corrections), the parton densities in a nucleus are quite generally given by 

F i/A (x, M)=Z F i/P (x, M) Sf) p (x, M) + (A- Z) F i/n (x, M) Sf) n (x, M), (2.5) 

in the perturbative calculations. Here, Sfj p (Sf^ n ) denotes the ratio of the proton 
(neutron) PDF in a nucleus A over that in a "free" proton (neutron). It is given 
throughout this paper by a global fit of Drell-Yan and DIS data on nuclear targets 
(EKS, [18]). Note that at very small x, however, coherence effects become large and 
the above factorised form for the nuclear PDFs may no longer be valid. Nevertheless, 
the typical x values we shall meet in this study are pretty large, say x > 0.01, hence 



the use of Eq. (|2.5|) appears to be fully justified. 

Finally, the partonic energy loss process in the dense medium is addressed. We 
shall follow a standard procedure, and model the energy loss through a shift in the 
z variable entering the fragmentation functions, [19] 

z^D™ d /k (z^,M F ,kJ = f k±il *"* de V k (e,kJ z* D^ /k (z%, M F ), (2.6) 

Jo 

Vd{t,k ± ) standing for the probability that the hard parton k with transverse mo- 
mentum k ± has lost an energy e. It was introduced in [20] and computed numeri- 
cally in [21] from the Baier-Dokshitzer-Mueller-Peigne-Schiff medium- induced gluon 
spectrum [22]. This spectrum - and thus eventually the quenching factors - de- 
pends on one scale, u> c . Accounting for the expansion of the produced medium, it 
reads [23,24] 

u c ~ q(t ) t L, (2.7) 

assuming a purely longitudinal expansion. The initial-time transport coefficient q(t ) 
in Eq. (|2.7[) represents the scattering power of the medium, and L the length covered 
by the hard parton in the medium. The transport coefficient has been determined 
perturbatively in a cold nucleus, q ~ 0.045 GeV 2 /fm [22], in good agreement with a 
phenomenological estimate from Drell-Yan production off nuclei [25]. In hot quark- 
gluon plasma, it is expected to scale like [24] 

q ~ a s e 3 / 4 , (2.8) 

where e is the medium energy- density and a s the strong coupling constant evaluated 
at a scale O (qL) . The energy- density reached in central heavy- ion collisions, esti- 
mated in the saturation model of Ref. [26], shows a parametric dependence e ~ Q 4 , at 
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an initial time to ~ 1/Q S , with Q s the gluon saturation-momentum. From Eqs. (|2.7| ) 
and ( [2~8l) , u c should scale like Q 2 S . From the v /s NN dependence of the saturation 
momentum [26], we thus have 3 



neglecting the running of a s in Eq. (|2.8|) . We shall take the uj c = 10-15 GeV range to 
predict the pion and photon quenching at ^/s NN = 62.4 GeV, to be consistent with 
our previous estimate, u c = 20-25 GeV, at RHIC full energy [8]. 

To summarise, quenching factors are computed, in the following, under various 
assumptions concerning the nuclear effects: 

(i) Only isospin corrections, Eq. (|2.4j ), are considered, that is, neglecting any nu- 
clear corrections to the PDFs nor energy loss effects. This prediction should 
be rather model- independent, and is labelled "Au Au" in the figures; 

(ii) Isospin and nuclear PDF (EKS) corrections are included in the nuclear parton 
densities, Eq. (|2.5| ). These predictions are labelled 'Au Au EKS"; 

(iii) Parton densities include isospin/nPDF corrections, Eq. ( |2.5| ), and medium- 
modified fragmentation functions, Eq. (|2.6|), are used to account for energy 
loss processes (labelled 'Au Au EKS u c = 10-15 GeV"). 

Although the above-mentioned mechanisms are expected to be dominant, it is 
worth to mention that other nuclear effects could affect the quenching of high-p x 
particles. It is the case for instance of the initial-state multiple scattering of partons 
in the nuclei, leading to the well-known Cronin effect (see e.g. [27] for a review). It 
is is expected to be significant at not too large p ± (say, below 5 GeV [28]) and is 
neglected in this study. 

3. Results 

3.1 Absolute spectra 

The transverse momentum spectra in p-p collisions at y/s NN = 62.4 GeV are shown 
in Fig. |l] for pions (left) and photons (right). The fragmentation process in the 
pion channel makes the cross sections decreasing much faster with p ± than what is 
predicted for the photons, which are produced directly in the partonic process at 
large x x = 2p ± /^/s NN > 0.1. Both processes turn out to have similar cross sections 
around p x ~ 15 GeV. 

3 Also interesting is the atomic mass dependence of the u> c scale. From [26], we get lo c ~ A 0/25 L ~ 
A 58 . This could be useful to predict the pion/photon quenching in lighter systems such as Cu-Cu 
at RHIC or Ar-Ar at the LHC. 




(2.9) 
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The theoretical uncertainty in the pion sector is pretty large, from a factor 2 
at low p ± up to a factor 3 at p x ~ 20 GeV. At large p ± , most of this uncertainty 
can be attributed to the fragmentation process (hence to the fragmentation-scale 
dependence) since fragmentation functions are probed at large z. On the contrary, 
it is interesting to remark that large-p ± photons, being produced directly, show a 
somewhat lesser scale-dependence (factor of 2 uncertainty over the whole p x range). 
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Figure 1: Pion (left) and photon (right) invariant cross section in p-p collisions at 
yfs = 62.4 GeV computed at NLO accuracy, varying simultaneously the factorisation, 
the renormalisation and the fragmentation scales from p ± / \/8 to p ± . Pions and photons 
are produced in the [—0.35; 0.35] rapidity interval. 



Pions and photons have been measured a long time ago in p-p collisions at y 7 ? = 
62-63 GeV at the ISR (see [29] for a list of references), and subsequent perturbative 
calculations were shown to describe fairly these data. We refer in particular the 
reader to Refs. [6,11,12] for a complete review of single-inclusive pion and photon 
QCD predictions from fixed-targets to collider experiments. Recently, all ISR 7r°-data 
were fitted with an empirical law [29] reproduced as a dashed line in Fig. |T[ The 
agreement between this fit and the present calculation is particularly good, although 
it is worth to remark that the "data" are best reproduced 4 with the present lower 
estimate, i.e. using all scales equal to p ± . 

3.2 Quenching factors 

Quenching factors in 0-5% central Au-Au collisions at y/s NN = 62.4 GeV are com- 
puted in Fig. [| for pions and photons. Let us first comment on the pion case (Fig. ^|, 
left). Isospin corrections prove tiny, with a 5% enhancement {R AuAu > 1), at most, 

4 It should however be kept in mind that the ISR measurements have a rough 20%-spread around 
this best fit. 
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at high p ± ~ 15 GeV. More pronounced are the nuclear effects in the nPDFs which 
deplete somehow (20%) the pion yield in nuclei due to the EMC effect at large 
x > 0.2. The slight enhancement seen at p ± ~ 5 GeV is due to the rather strong 
gluon anti-shadowing assumed in the EKS parametrisation around x ~ 0.1. Overall, 
both isospin and nPDFs corrections remain small and do not affect the pion quench- 
ing factor by more than 20% over the whole p ± range. These prove, however, more 
pronounced at ^/s NN = 62.4 GeV than at 200 GeV due to the larger x probed in the 
nuclear parton densities. 




o 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 o 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 

5 10 15 5 10 15 

Pt (GeV) p T (GeV) 



Figure 2: Pion (left) and photon (right) quenching factors in central Au-Au collisions at 
^/s NN = 62.4 GeV. Calculations are done at LO, assuming (i) isospin (solid), (ii) isospin and 
nPDFs (dotted), (iii) isospin, nPDFs and energy loss (band) effects (see text for details). 

On the contrary, the energy loss mechanism quenches the pion p ± spectrum by 
roughly a factor of 5! Furthermore, the quenching factor is seen to be rather p x - 
independent; this flat behaviour coming from the interplay of nPDFs corrections at 
large x (EMC effect) and energy loss processes. Quite surprisingly, the quenching 
turns out to be as pronounced as at y/s NN = 200 GeV, despite the twice smaller up- 
scale. The reason is actually due to kinematics. With decreasing ^/s NN and at a given 
p ± (i.e. with increasing x ± ), the phase-space available to produce high-/c ± partons 
becomes more and more restricted. The fragmentation process into the pion hence 
occurs at larger z, where FFs are known to fall dramatically, and therefore where 
energy loss effects should be the strongest. This - as well as the more pronounced 
nPDF corrections - explains why the pion quenching factor proves as small as ~ 0.2 
despite the much smaller uj c = 10-15 GeV assumed in the calculation. 

The pion quenching at y/s NN = 62.4 GeV has been already computed within 
various parton energy-loss models by Vitev [28], Adil and Gyulassy [30], Wang [31], 
and Eskola, Honkanen, Salgado, Wiedemann [32]. All calculations predict slightly 
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larger quenching factors, FT ~ 0.3-0.4, in better agreement with PHENIX prelim- 
inary results in Au-Au collisions [33] normalised to ISR p-p data [29]. It should 
be interesting to see whether the final PHENIX p-p measurements confirm this first 
observation. The overestimate of the quenching in the present model at low p ± ~ 4- 
6 GeV has already been pointed out at y/s NN = 200 GeV in Ref. [8] and is possibly 
due to the Cronin effect not considered here 5 . Hopefully future data will be able to 
tell whether this overestimate subsists at larger p ± . 
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Figure 3: Ratio of the overall photon production over the background, 7 all /7 bkgd , in central 
Au-Au collisions normalised to the p-p case. LO calculations are done in p-p and Au-Au 
collisions at ^/s NN = 62.4 GeV. 

The quenching pattern in the photon channel (Fig. ^ right) is quite different. 
Already noticeable in the predictions at y/s NN = 200 GeV [8], and possibly seen in the 
PHENIX preliminary measurements [7], the isospin corrections at y/s NN = 62.4 GeV 
now prove spectacular: from p ± = 4 GeV to p ± = 15 GeV, the quenching factor 
decreases from 0.85 down to 0.65. Because of the electromagnetic coupling, the 
scattering of up-quarks is indeed favoured as compared the down-quarks due to their 
respective electric charge. The relative lack of valence up-quarks in the neutron 
(as compared to a proton) tends to deplete photon cross section in large nuclei 
(with Z ~ A/2), which explains why the quenching factor is expected to be smaller 
than 1. Nuclear PDF corrections further increase the negative slope of the photon 

5 We may also remark that the energy dependence assumed here, Eq. Q2.9Q , is slower than the 
y/s^ ' 57 assumed in Ref. [32]. 
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quenching, making the ratio R ~ 1 at lowp x because of anti-shadowing, and R ~ 0.5 
in the highest p ± bin (EMC effect). As emphasised in [8], energy loss could also 
affect prompt photon production in Au-Au collisions because of the fragmentation 
contribution (see Eq. ( |2.2|) ). However, its effects remain rather small, roughly 20% 
for all p ± . As a matter of fact, at low y/s NN , photons are mostly produced directly 
in the hard subprocess and should not be too much quenched by the medium. In 
that sense, the restricted phase-space has opposite consequences - as far as energy 
loss processes are concerned - in the pion and the photon channel, because of their 
different dynamical production process. In particular, neglecting any isospin/nPDF 
corrections, the quenching factor at the very edge of phase-space should be, strictly 
speaking, R n (2p ± = y/s NN ) = for pions but R n (2p ± = y/s NN ) = 1 for photons. 

This difference is well illustrated in the 7 all /7 bkgd ratio, that is the overall photon 
yield (including the prompt photon signal and the n° — > 77 decay background 6 ) 
normalised to the background. This ratio is, however, quite sensitive to the absolute 
magnitude of the 7r° and 7 spectra in p-p collisions. In order to lower this uncertainty, 
Fig. I displays the double ratio (7 all /7 bkgd ) AuAu / (7 all /7 bkgd ) pp as a function of the 
transverse momentum. This observable looks promising as it is clearly sensitive 
to energy loss processes - see e.g. the strong enhancement in Fig. [| and its u> c 
dependence - while rather independent of isospin/shadowing corrections when going 
from p-p to Au-Au collisions. 

3.3 Comparison with y/s NN = 200 GeV 

Naively, we could think that the smaller the incident energy y/s NN (and therefore the 
smaller the medium energy-density), the weaker the quenching. Looking at Eq. ( |2.9| ). 
this would mean that the quenching factor smoothly decreases with ^/s NN . This is 
not necessarily true, for two reasons. The first one comes from the nPDF corrections: 
parton density ratios, S^, clearly have a non-monotonic behaviour as a function of 
Bjorken-x or, alternatively, as a function of y/s NN at a given p ± . The second reason 
deals with the phase-space restriction at small y/s NN > u c , already mentioned in 
Sect. |3]^: energy loss effects can be magnified at low y/s NN because of the larger z 
values probed in the fragmentation process (or equivalently because of the steeper 
partonic spectra at large k ± < 

In order to compare the quenching at intermediate and full RHIC energy, the 
double ratio 

i? AA Mw = 62.4 GeV) 

K 62.4/200 — -5 / . onn n \r\ ' ^ 6 

Raa{t/s^ = 200 GeV) 

is computed in Fig. £|. This ratio has been introduced by Adil and Gyulassy in [30] 
in order to get rid of the uncertainty in the initial gluon rapidity-distribution, and 
was determined in the pion channel assuming energy loss effects only. 

6 As in Ref. [8], the tt° decay spectrum is estimated to be 2/(n — 1) that of the single-inclusive 
7T° spectrum, where n represents the power-law exponent of the 7r° spectrum in p-p collisions [4] . 



- 9 - 



Although the energy loss probability distribution used here is somehow differ- 
ent to the one assumed in Ref. [30], it is worth to mention that roughly the same 
trend is observed in our model (Fig. f|, left), and especially the negative p ± -slope of 
this ratio coming from the phase-space restriction above discussed. Our prediction 
is in particular similar to their calculation assuming an initial gluon distribution 
dN g /dy = 770 [30]. The inclusion of nPDF corrections only affects the magnitude of 
the ratio, lowering it by 20% (dashed-line). Consequently, the transverse momentum 
above which the quenching is larger at v /s NN = 62.4 GeV than at y/s NN = 200 GeV 
(R 6 2.4/2oo(pJ < 1) is p ± ~ 8 GeV instead of p ± > 12 GeV in Ref. [30]. 




P T (GeV) p T (GeV) 

Figure 4: Ratio of the quenching factor at ^/s NN = 62.4 GeV over that at ^/s NN = 200 GeV 
in the pion (left) and the photon (right) channel. Calculations are done at LO, assuming 
(i) isospin (solid), (ii) isospin and nPDFs (dash-dotted), (iii) isospin, nPDFs and energy 
loss (band) effects (see text for details). 

The same calculation is carried out in Fig. |4] (right) in the photon channel. 
Despite the different dynamical production process, the ratio R 62.4/200 including 
energy loss is remarkably similar to what is found in the pion sector. This is, however, 
largely accidental, since the isospin corrections in the double ratio turn out to be 
really different for pions and for photons. Interestingly, the double ratio is sensitive 
to energy loss effects at low p ± < 7 GeV. At larger p ± > 11 GeV, the ratio does 
not depend much on the energy loss mechanism any longer (because of the direct 
process) but rather should probe the EMC effect in the nPDFs. 

4. Summary 

Large-p x pion and photon spectra have been considered in p-p and Au-Au scattering 
at intermediate RHIC energy, y/s NN = 62.4 GeV. After discussing NLO predictions 
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in p-p collisions, quenching factors in the 0-5% most central Au-Au collisions are 
predicted for both pions and photons within the same model. The energy loss process 
quenches rather strongly the pion yield, but proves rather inefficient for the photons, 
whose quenching factor is dominated by significant isospin and nPDF corrections. 
The quenching of the photon total yield over its background is also determined and 
turns out to be rather sensitive to the energy loss processes. 

Predictions at y/s NN = 62.4 and at y/s NN = 200 GeV are then compared. Despite 
the smaller gluon density, the large-p ± pion quenching at ^/s NN = 62.4 GeV proves 
actually as pronounced as at full RHIC energy because of the strong phase-space 
restriction at this energy, as already mentioned in [30,32]. On the contrary, the 
kinematic dependence of the photon quenching is opposite because of the direct 
contribution which dominates the inclusive yield: when the phase-space becomes too 
restricted - i.e. at low y/s NN -, photons are produced directly and should therefore 
not be affected by the dense QCD medium. 
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